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Abstract. In this article we present numerical results of conduction in a disordered 
quasi-lD wire in the possible presence of magnetic impurities. Our analysis leads us to 
the study of universal properties in different conduction regimes such as the localized 
and metallic ones. In particular, we analyse the cross-over between universality classes 
occuring when the strength of magnetic disorder is increased. For this purpose, we 
use a numerical Landauer approach, and derive the scattering matrix of the wire from 
electron's Green's function 
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1. Introduction 

Interplay between disorder and quantum interferences leads to one of the most 
remarkable phenomenon in condensed matter : the Anderson localization of waves. 
The possibility to probe directly the properties of this localization with cold atoms [1, 2] 
have greatly renewed the interest on this fascinating physics. In this paper, we focus on 
the particular situation where electrons encounter two kinds of disorder: a usual scalar 
potential at the origin of diffusion, and a magnetic potential, arising from a collection of 
frozen random magnetic moments. This situation is naturally realized experimentally 
in the study of transport properties of metallic spin glass wires [3, 4, 5, 6]. In these 
wires, the spins freeze at low temperatures when entering the spin glass phase due to 
the frustrating magnetic couplings. In this glassy phase, and neglecting any residual 
Kondo effect in this regime, the impurities act effectively as a (weak) magnetic potential. 
We study numerically the effect of both types of disorder on the statistical properties 
of the wire conductance. In particular, we will focus on the experimentally relevant 
crossover of (weak) localization properties of the wire as a function of the magnetic 
disorder strength. 

One dimensional disordered electronic systems are always localized. Following 
the scaling theory [7] this implies that by increasing the length L x of the wire for a 
fixed amplitude of disorder, its typical conductance ultimately reaches vanishingly small 
values. The localization length £ separates metallic regime for small length L x <C £ from 
the asymptotic insulating regime. In the present paper, we focus on several universal 
properties of both metallic and insulating regime of these wires in the simultaneous 
presence of two kinds of disorder. The first type corresponds to scalar potentials induced 
by the impurities, for which the system has time reversal symmetry (TRS) and spin 
rotation degeneracy. In this class the Hamiltonian belongs to the so-called Gaussian 
Orthogonal Ensemble (GOE) of the Random Matrix Theory classification [8] (RMT), 
corresponding to the class AI in the modern classification of Anderson universality 
classes (see e.g. [9]). If impurities do have a spin, the TRS is broken as well as spin 
rotation invariance. The Hamiltonian is then a unitary matrix, which corresponds 
in RMT to the Gaussian Unitary Ensemble (GUE) with the breaking of Kramers 
degeneracy [10], and to the Anderson class A [9]. However, for the experimentally 
relevant case of a magnetic potential weaker than the scalar potential, the system is 
neither described by the GUE class, nor by the GOE class, but extrapolates in between. 
This intermediate regime, of particular relevance experimentally, is the main object 
of study of the present paper. Moreover the present work paved the way towards a 
numerical study of the correlation of conductances in the cross over regime [11]. 
This paper is organized as follows : in section 2, the model and the numerical method 
used will be described in details. In section 3 we identify the localized and metallic 
regime of transport of the system. The localization length will be determined by two 
different methods and the cross over between both universality classes (GOE, GUE) 
will be highlighted. In section 4, the insulating regime is studied with a particular focus 
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on the statistical distribution of the conductance, which allows us to highlight universal 
behavior. In section 5, we turn to the study of the metallic regime, and perform a 
careful analysis of each of the first three cumulants of the statistical distribution of the 
conductance. We focus on the universal properties of conductance fluctuations, and the 
non-analyticity of the complete distribution is discussed. Finally section 6 is devoted to 
the conclusion. 

2. The model and the method 

2.1. The model 

In this paper, we study numerically the scaling of transport properties of wires in 
the presence of magnetic and scalar disorders. We will focus on the regime of phase 
coherent transport, reached experimentally at low temperature (see in particular [6]). 
In this regime, the phase coherence length which phenomenologically accounts for 
inelastic scattering of electrons on impurities [12] is larger than (or comparable with) 
the wire's length L x , so that phase coherence for the propagating electrons is preserved 
in the whole sample. Note that this phase coherence length includes in particular 
a contribution from inelastic scattering on the non frozen magnetic impurities through 
a Kondo dephasing, which is strongly reduced when entering the magnetic glass phase 



We describe the behavior of electrons inside the disordered wire using a tight- 
binding Anderson lattice model with two kinds of disorder potentials: 



tij is the hopping term from site i to j. In the following, tij will take two different values: 
tij = t/j in the longitudinal x direction and t^ = t± in the transverse y direction. The 
scalar disorder potential V = {v{\i is diagonal in electron-spin space. We choose the 
Vi to be random scalars uniformly distributed in the interval [—W/2,W/2]. In this 
work, we have chosen without loss of generality to fix t// = 1 so that all energy scales 
are relative to the bandwidth t// = 1, and the amplitude of disorder W = 0.6. In 
eq. (1) s, s' label the SU(2) spin of electrons and the Si account for spins of the frozen 
magnetic impurities. A realistic choice for these frozen spins in a spin glass phase 
consist in considering classical spins with random orientations, thereby neglecting any 
small spatial correlation in a spin glass configuration [13]. The coupling J between the 
electron spins and the magnetic impurities fixes the amplitude of the magnetic disorder. 
In this work, we will monitor the behavior of the transport properties of the sample as a 
function of this amplitude J, from J = to J = 0.4. Indeed, variation of the amplitude 
of magnetic disorder J allows to extrapolate from GOE / class AI for J = to GUE / 
class A for J ^ 0. We will also show in section 5 that, as a bonus, the presence of the 
magnetic disorder allows also an numerically easier settlement of the universal metallic 
regime of the wire. 



[6]. 




(1) 
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For a given realization of both scalar {vi}i and magnetic disorder {Si}i, the 
Landauer conductance of this model on a 2D square lattice of size L x x L y is evaluated 
numerically using a recursive Landaueur method described in details in the next 
paragraph. 

2.2. The method 

We have chosen to use a numerical method based on the lattice model (1) as opposed 
to e.g. random matrix or Dorokhov-Mello-Pereyra-Kumar (DMPK) [14, 15] method so 
as to provide a numerical study allowing for an easy comparison with the experimental 
situation[6]. Moreover, this method allows for further natural developments such as the 
study of the conductance change upon magnetic impurities spin flipping, which would 
be difficult to reach by alternative method. Starting from a lattice model such as (1), 
the natural method providing the conductance of a finite size sample is based on the 
Landauer formalism [16]. 

We consider a two-terminal setup, with electrodes connected to the wire at x = 
and at x = L x . These electrodes are described as semi-infinite ribbons with the same 
transverse geometry as the sample, and described with (1) but without randomness. 
Electrons are then confined in the transverse y-direction via a potential that has the 
form: 

V(y) = if < y < L y ; V(y) = oo otherwise (2) 

This confining potential in the y direction leads to the appearance in the electrodes of 
N modes propagating in the x-direction. The complete wave function of an electron in 
this tight binding lattice model reads then: 

^(x,y) = My)e tk * x , (3) 
where k x is the momentum of electrons in the longitudinal direction and 



We used N y = L y in units of lattice spacing. The group velocity of this mode reads: 

v n = 2^ sin (k x ) . (5) 

where t// is the longitudinal hopping amplitude. This velocity depends on the 
momentum k x which is determined for a constant energy by the dispersion relation 

(T17T \ 
N - 1 J - 2t// cos(k x ). (6) 

At given energy E — (j,, we end up with the following relation for the longitudinal part 

of the momentum of the electron: 

(H-E t ± ( nn \\ 
k x = arccos — cos — . (7) 



Cross- Over between universality classes in a magnetically disordered metallic wire. 5 



To optimize the efficiency of the numerical study, we fix t// = 2t± and stay near the 
band center, avoiding in particular the presence of fluctuating states studied in [17]. 

To compute the conductance of such a wire, the Landauer-Biittiker formalism of 
coherent transport is used [18]. It allows to relate the dimensionless conductance g of a 
diffusive wire with the scattering matrix T: 

9= Tmn ' ( 8 ) 

modes m,n 

where T mn is the transmission coefficient between modes to and n. The dimensionless 
conductance g is defined from the conductance G of the system as g — Gj (2e 2 /h) when 
spin degeneracy is present (J = 0), and g = G/(e 2 /h) otherwise (J ^ 0). 

Following Fisher and Lee [19] we relate the scattering matrix to the electronic 
retarded Green's functions of the system through: 

Ny 

tmn = lhy/V n V m ^ ( l ) n(y)Q R (y, X = (%', X = L x )(f) m (y'). (9) 
J/,2/'=l 

with T = Tr(tH), v n and n (resp. v m and <f> m ) are the group velocity and the eigen 
wave function of propagating mode n (resp. to). Mode n belongs to the left lead 
whereas mode to belongs to the right one. In (9) Q R (y,x = 0\y',x = L x ) represents 
the retarded Green's function of an electron between the point (x = 0,y) in the left 
contact between the system and the electrode and the point (x = L x ,y') in the right 
contact. Consequently, by calculating the retarded Green's functions of the system only 
between its both sides, we are able to determine the dimensionless conductance of the 
wire. An efficient method to calculate Q R (y, x = 0\y r , x = L x ), which takes advantage of 
the quasi-one dimensional nature of the ribbon consists in obtaining it recursively [20]. 
In figure 1 the principle of the method is sketched: using a Dyson equation we deduce 
the boundary Green's function of a system of size n + 1 from the corresponding Green's 
function of a subsystem of size n, and the exact Green's function of the n + 1 row. This 
allows to perform matrix inversion only of the simple row system. At both initial and 
final steps, we reconnect the system to the semi-infinite electrodes (see fig. 2) described 
by a standard self energy: 

1 Xy 

Cund (2/1,2/2) = --J2Myiy kx Uy2). (10) 

The last step consists in combining the Green's functions of the sample of desired 
longitudinal size with the one of the right lead, which is also given by equation (10). 
The unitarity of the corresponding scattering matrix {T mn } m n is used to monitor the 
accuracy of the numerical method. Such test yielded typical relative error of order 10~ 4 
for a system size L x = 1600 and L y = 40. This method allows one to compute the 
conductance of a wire of length L x and of width L y for any given configuration of scalar 
disorder V and for any configuration of frozen classical spins 

In the next sections, we study the properties of this conductance for one given 
configuration of magnetic disorder but for many different configurations of scalar 
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G R (l,n) g(n+l,n+l) G R (l,n+l) 

Figure 1. Principle of recursive calculation of retarded Green's functions of the wire. 
Use of a Dyson: G^ +1 = G,f + G^VG^ +1 . At each step the longitudinal length L x is 
increased by one lattice spacing. 



Figure 2. Boundary conditions: the wire is connected to two leads represented by 
two semi infinite metallic wires. 



disorder. Universal properties are identified by varying the transverse length L y from 
10 to 80, with the aspect ratio L x /L y taken from 1 to 6000. Typical number of 
configurations of scalar disorder V we used were = 5000, with exceptions for the 
study of the localization properties where for L y = 10 we sampled the conductance 
distribution for 50000 different configurations of disorder. 

3. Localization length 

3.1. Determination of £ 

We start our analysis by a determination of the parameters corresponding to the metallic 
and localized regimes, through a careful determination of the localization length of the 
system. While experimentally the only accessible regime of a phase coherent wire is the 
metallic regime, numerically this regime is difficult to reach and describe quantitatively, 
as opposed to the localized regime. This is related to the extreme reduction in the 
number of propagating modes in the numerical system which is associated with a 
corresponding reduction of the localization length. Hence in order to clearly identify the 
conditions to access universal properties of the transport in the metallic regime, we start 
by a detailed determination of this localization length in the experimentally relevant 
crossover situation. Afterwards we will take the opportunity of the present study to 
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describe other characteristics of the localized regimes in the crossover situation, before 
turning to our main interest : the universal metallic regime. 

The localization length separates short wires of metallic behavior from a insulating 
long wires. A first method to access to the localization length £ from the conductance 
consists in considering the scaling behavior of the typical conductance g typ defined as: 

9t yp = e (log9> , (11) 

where (•) represents the average over the different configurations of scalar disorder 
V. This typical conductance decays exponentially with the longitudinal length of the 
wire [21, 22] : 

2L X 

9ty P ~ e e (12) 
in the regime of long wires L x 3> £. Figures 3, 4, 5 and 6 show the behavior of 




x 10 



Figure 3. Evolution of (log g) as a function of longitudinal size for different transverse 
length (L y = 10,20,30,40) and J = 0. The linear part of the curve allows one to 
get the localization length £(J = 0) from the scaling form in the insulating regime 



(log g) as a function of longitudinal length L x for different widths L y . Different curves 
correspond to different values of magnetic disorder. The linear fit of the large length 
part of the curves allow for a precise determination of the corresponding localization 
length for each value of L y and J. 

A second method to determine this localization length consists in considering the 
Lyapunov exponent 7 of the transfer matrix of the system, following a standard random 
matrix theory approach[23, 21, 24]. This exponent can be deduced from the conductance 
as 



1 

2T 



log 1 



(13) 
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Figure 4. Evolution of (log g) as a function of longitudinal size for different transverse 
length (L y = 10,20,30,40) and J = 0.05. The linear part of the curve allows one to 
get the localization length £(J = 0.05) from the scaling form in the insulating regime 
(logff) 
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Figure 5. Evolution of (log g) as a function of longitudinal size for different transverse 
length (L y — 10,20,30,40) and J = 0.2. The linear part of the curve allows one to 
get the localization length £( J = 0.2) from the scaling form in the insulating regime 
(log S > = -*4*. 



and the localization length follows from its asymptotic behavior : 

r 1 = Hm 7 (L X ). (14) 

L x — ¥oo 

On figure 7 and 8 we have plotted the Lyapunov exponent versus the inverse of the 
longitudinal length for different values of magnetic disorder. Different curves correspond 
to different widths of the wire. With this method, a simple extrapolation of the curve 
is necessary to obtain £, without any fit. the value of the inverse of the localization 
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Figure 6. Evolution of (log g) as a function of longitudinal size for different transverse 
length (L y — 10,20,30,40) and J = 0.4. The linear part of the curve allows one to 
get the localization length £( J = 0.4) from the scaling form in the insulating regime 
(logff) 
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Figure 7. Evolution of Lyapunov exponent with the inverse of longitudinal length 
in semi-log plot. Circles correspond to J = 0, squares to J — 0.05, diamonds to 
J = 0.2 and triangles to J = 0.4. The value of the transverse length is L y = 10. The 
localization length can be extrapolated from the value of 7 for L x — > 00. 



length for an infinite wire. Nevertheless, we have found that this method shows less 
accuracy than the preceding one: as shown on figures 7 or 8, the Lyapunov exponent is 
still varying for the longest longitudinal length. We find that both methods give fully 
compatible results while the Lyapunov exponent method requires much larger system 
sizes than the typical conductance method for a given required accuracy. 

A first manifestation of the universality of the Anderson localization classes appears 
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Figure 8. Evolution of Lyapunov exponent with the inverse of longitudinal length 
in semi-log plot. Circles correspond to J = 0, squares to J — 0.05, diamonds to 
J = 0.2 and triangles to J = 0.4. The value of the transverse length is L y = 20. The 
localization length can be extrapolated from the value of 7 for L x — > 00. 



through the dependance of £ on the transverse length L y (or the number of propagating 
modes). It is expected to follow[21, 25]: 

i = ((3L y + 2 - (3%, (15) 

with l e the elastic mean free path and (3 encodes the universal class of the model : 
(3 = 1 corresponds to the orthogonal universality class GOE while (3 = 2 for GUE. 
Note that this change in [3 is accompanied by an artificial doubling of the number of 
transverse modes N y = L y — )■ 2N y due to the breaking of Kramers degeneracy [21]. 
This effective factor 4 when breaking the spin rotation symmetry has been discussed in 
ref. [26, 27] in details, when discussing the magnetic field dependance of this localization 
length, in comparison with random matrix and Non-linear sigma models. Comparison 
of numerical localization lengths for different J with (15) is shown in fig 9. Excellent 
agreement is found for J = (GOE class, (3 = 1) and with the GUE class for J > 0.2. 
For intermediate values of J 7^ we observe a crossover between the two extreme GOE 
and GUE laws, which cannot be described by eq. (15), and for which no analytical work 
exists to our knowledge. 

From these results, we also notice that the localization regime is reached for much 
longer wires in the presence of magnetic impurities (GUE case) than without (GOE): 
localization is hampered by the presence of these magnetic impurities. As we will discuss 
below, this property helps in observing numerically the universal weak localization 
regime and the associated universal conductance fluctuations. 
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Figure 9. Evolution of localization length as a function of transverse length. l e is the 
mean free path of the diffusive sample. Different behavior of the localization length if 
J = or J ^ 0. Inset : Scaling of the typical conductance (log g) 
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3.2. The Insulating and Metallic regimes 

The localization length discriminates between both insulating and metallic regimes: the 
ribbon behaves indeed as a metal (g 3> 1) for lengths L x <C £ and as an insulator 
(g <C 1) if L x 3> £. In both asymptotic regimes the shape of the Probability Density 
Function (PDF) of the conductance which is known: it is Log-normal for insulating 
wires[21] and gaussian for metallic ones. By varying the longitudinal length we can 
study the evolution of this PDF from a gaussian to a log-normal distribution, as shown 
on figure 10. This plot is done for a given value of L y and J. One can notice that 
in the metallic regime, the PDF is very well approximated by a gaussian for relatively 
small wires : the gaussian regimes is easily reached, whereas it takes length much larger 
than the localization length for the distribution to become log-normal in the insulating 
regime. This point will be discussed more precisely below on the cumulants of this PDF. 
The insulating regime is then characterized by (g) < 1 and the metallic one by (g) > 1. 

In order to characterize samples by the average (g), and in particular plots higher 
cumulants as a function of (g), we now turn to a short study of the behavior of this 
first cumulant as a function of the system size. On figure 11 we have plotted (g)(L x ) 
for different values of magnetic disorder J (hence different localization lengths). These 
curves approximately collapse when plotted against the scaling variable L x /£, as show 
on fig. 12. we remind the reader that for J = 0, the average conductance is defined by 
g = G/2Gq, which explains why J = and J = 0.4 curves coincide on fig. 11: according 
to figure 9, only for these values of magnetic disorder universality classes are reached. 

This study of (g)(L x ) allows to proceed in the study of higher cumulants of the 
PDF of g and test prediction of the single parameter scaling of distributions [28]. 



Cross- Over between universality classes in a magnetically disordered metallic wire. 12 




J=0, L = 10, N = 5 0000 points 

' _Num data, L 

X 

_Num data, L 

X 

Num data, L 

_ X 

_Num data, L 



400 
800 

: 1600 

2000 



— Gaussian interpolation 
---Gaussian interpolation 
• Log-normal interpolation 
— Log-normal interpolation 



1-5 g 



Figure 10. Evolution of the statistical distribution of the conductance for different 
longitudinal sizes for L y = 10 and J = 0. N is the number of disorder configurations 
used. Plain lines are gaussian (if (g) > 1) or log- normal (if (g) < 1) interpolations of 
numerical data. 
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Figure 11. Evolution of average conductance versus longitudinal length for different 
values of magnetic disorder for L y = 10 and J = 0,0.05,0.1,0.2 and 0.4. The line 
(g) = 1 is plotted as a frontier between insulating and metallic regimes. 



4. Universal Insulating Regime 

4-1. Probability Density Functions 

In the insulating regime L x > £, we expect a Log-normal conductance statistical 
distribution, as seen previously. However in the weakly insulating regime (g) < 1 
this log-normal asymptotic form is not reached. Instead, as shown in fig. 10 we find a 
non-analytical behavior of P(g) in agreement with [29, 30, 31, 32]. In order to study 




Figure 12. Evolution of average conductance versus the scaling variable L x /£, for 
L y = 10 and J = 0, 0.05, 0.1, 0.2 and 0.4. All curves nearly collapse in one single curve. 
The frontier between insulating and metallic is drawn again. 




Figure 13. Comparison of Probability density functions (PDF) of conductance for 
,/ = (plain curves) and J = 0.2 (dashed curves). Plots are performed for different 
values of average conductance, (a): (g)(J = 0) = 0.84 and (g)(J = 0.2) = 0.79. 
(b): (g)(J = 0) = 0.67 and (g)(J = 0.2) = 0.62. (c): (g)(J = 0) = 0.45 and 
(g)(J = 0.2) = 0.42. (d): (g)(J = 0) = 0.21 and (g)(J = 0.2) = 0.18. 



the dependance of this non-analyticity on the universality class, we have plotted on 
figure 13 the distribution P(g) for similar values of (g) but different magnetic strengths 
J = (GOE) and J = 0.2 (GUE). The shapes of these distributions are highly similar 
if (g) <C 1, showing that distributions for J = and J ^ reaches the same Log- normal 
distribution at large system sizes, in agreement with the super-universality scenario[33]. 
In the intermediate regime ((g) ~ 1), shapes are symmetry dependent. Moreover we find 
that the non-analyticity appears for different values of conductance (close to 1) and the 
rate of the exponential decay [29] in the metallic regime seems to differ from one class to 
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the other (see for instance curves (a) or (b)). Finally, the plot on figure 14 represents the 
distribution of conductance of mean value just above and below the threshold (g) — 1. 
Plain lines represent gaussian interpolations with a mean and a variance given by the 
first and the second cumulant of each numerical conductance distribution. For (g) > 1, 
the gaussian interpolation approximates very well the full distribution while as soon as 
(g) < 1, the gaussian approximation only applies in the tail g > 1 of the distribution of 
conductance [32]. The shape of the distribution for g < 1 (Figure 13) agrees qualitatively 
with the numerical results of [34] (see in particular their fig. 4), Unfortunately, a more 
accurate comparison proves to be difficult due to the lack of analytical description of 
the distribution. To quantize further these results on the whole distribution of log g we 
now turn to a quantitative study of the second and third cumulant of this distribution. 

4-2. Study of cumulants 

This conductance distribution converges to the Log-normal only deeply in the insulating 
regime, the convergence being very slow (much slower than in the metallic regime). This 
qualitative result is confirmed by the study of moments: in the insulating regime the 
second cumulant is expected to follow[35, 25]: 

((logs - (logs)) 2 ) = ((logs) 2 ) c = -2(log<?>, (16) 

Our numerical results are in agreement with this scaling with however very slow 
convergence towards this law: corrections are measurable even for the largest system 
size where the system is deeply in the localized state, as shown on fig. 15. More precisely, 
we find that for the deep insulating regime ((logg) 2 ) c = — 1.88(logg) slope —1.88, with 
a slight discrepancy with (16). This plot also shows that in the deep insulating regime 
the behavior of the second cumulant as a function of the first one does not depend on 
the value of magnetic disorder: both curves for J = and J = 0.2 follow the same 
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Figure 15. Plot of the variance of logg as a function of the mean for the orthogonal 
(J = 0) and unitary (J = 0.2) case, for L y = 10 and 20. The slope of the linear fit is 
equal to —1.88. This plot shows also super- universality as the behavior of the second 
cumulant does not depend neither on geometry of the wire nor on the universality 
class. 



law. This is in agreement with our previous result on statistical distributions: there is 
a super-universal behavior in the deep insulating regime. 

Finally we have studied the third cumulant of logg scaled in figure 16 as a function 
of the first cumulant. The linear behavior for each value of magnetic disorder in the 
deep insulating regime ((logg) < —4) is in agreement with the single parameter scaling. 
We find that contrary to the second cumulant the coefficient of proportionality between 
the skewness and the average depends on the symmetry of disorder, which denotes a 
lack of super- universality concerning this cumulant. For instance dots and diamonds 
(which correspond to the case J = 0) have the same behavior, as opposed to the case 
J = 0.2 (squares or triangles). A systematic study of this point with even larger system 
sizes and other numerical methods would be of high interest but is definitely beyond 
the scope of the present paper. 

The study of cumulants of the distribution of \ogg confirm the single parameter 
scaling of the distribution, with a slight discrepancy concerning the value of the 
coefficient of proportionality between second and first cumulant. Moreover, super- 
universality has been highlighted concerning the second cumulant but is lacking 
concerning the third one. 

4-3. Comparison with exact results in the cross-over regime 

For localization in wires connected to ideal contacts, exact formula for the average 
conductance [36] and conductance fluctuations [37] have been derived for the two 
universal orthogonal and unitary classes. These formula are of particular interest in the 
intermediate regime between the metallic (L x <C £) and the deeply localized (L x ^> £) 
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Figure 16. Plot of the skewness of log g as a function of the mean for the orthogonal 
(J = 0) and unitary (J = 0.2) case, for L y — 10 and 20. 



regime. We have numerically evaluated the formula (3.105) and (3.106) of ref. [37] and 
we compare them with our numerical data in Fig. 17 and Fig. 18. We find an excellent 
agreement between the J = data and the orthogonal exact formulae on one hand, and 
the J ^ data and the unitary formulae on the other hand. This comparison naturally 
breaks down for small sizes where the quasi-ld assumption for diffusion breaks down and 
a non universal regime takes place. As shown on Fig. 18, the exact unitary behavior is 
recovered for J = 0.05 beyond the cross-over length L m (see next section). Whether the 
scale dependance of (g) and (Sg 2 ) in the cross-over regime (J small such that L m ~ £) 
is amenable to an exact formula along the lines of ref. [37] is an open question of high 
interest. 



5. Universal Metallic Regime 

We now focus on the universal metallic regime described by weak localization. By 
definition weak localization corresponds to metallic diffusion, expected for lengths of 
wire l e <C L x <C £. For this regime to be reached, we thus need to increase £ through 
an increase of the number of transverse modes L y with all other parameters fixed (see 
eq. (15)). Moreover, for a fixed geometry, this regime will be easier to reach with 
magnetic impurities than without. As we saw on fig. 13, the shape of the PDF of 
conductance is a truncated gaussian in this regime. In the following we study its three 
first cumulants quantitatively, starting with the variance. 



5.1. Conductance Fluctuations and Universal Crossover 

In the weak localization regime, the conductance fluctuations (g 2 ) c are expected to be 
independent on the size of the system, and only depend on the universal localization 
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Figure 17. Comparison of the numerical average conductance (g) as a function of the 
reduced length L x /£ with the exact expressions of ref. [37] for the orthogonal (J = 0) 
and unitary (J = 0.2) case, for various transverse sizes L y . 
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Figure 18. Comparison of the numerical average conductance fluctuations (Sg 2 ) as 
a function of the reduced length L x /£ with the exact expressions of ref. [37] for the 
orthogonal (J = 0) and unitary (J = 0.4) case, for various transverse sizes L y . 

class of the model. The figure 19 shows that for a suitable value of transverse length, the 
system reaches the expect plateau in conductance fluctuations. The value of the plateau 
identifies with the expected values (1/15 and 4/15 for GUE and GOE respectively) with 
a high accuracy. However, the presence of this plateau depends strongly on the value 
of transverse length L y and on the magnetic disorder. Figures 20 and 21 illustrates this 
point with further details: these plots show the conductance fluctuations as a function 
of longitudinal size for J = and J = 0.2 and for two values of transverse length L y . On 
the first plot, for both values of magnetic disorder the universal plateau arises, whereas 
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it appears only for J = 0.2 if L y = 40. 

The evolution of this variance of the PDF of g depends on the longitudinal length 
through [12]: 

<tf> = {9% = if(0) + \f (*yf) + \f + If (*yf) ,(i7) 

where x = L x /L m and the scaling function F(x) depends only on dimension [12, 11]. 
The universality occurs in this equation in the two limit x = corresponding to J = 
or x 1 where F(x) — > 0. In both cases the variance becomes geometry independent. 
Moreover this expression shows that the whole crossover between the two classes is 
described by a universal crossover function parametrized solely by the length L m , called 
the magnetic dephasing length scale [12]. In the figure 22, these conductance fluctuations 
are plotted as a function of longitudinal length L x for different values of J. A single 
parameter fit by (17) provides the determination of the magnetic dephasing length L m 
as a function of the magnetic disorder J. The Universal behavior is also highlighted 
for strong enough magnetic disorder. The determination of this scattering length allows 
a precise study of average conductance, and in particular the study of the classical 
part as described in the next sub section below. Fig. 23 shows the scaling form of 
these fluctuations (as a function of L x /L m (J)) in excellent agreement with the theory 
(17). Moreover, for long wires (and large values of J) conductance fluctuations are no 
longer L x dependent and equal to 1/15, in units of Gq. This is the so-called Universal 
Conductance Fluctuations (UCF) regime which is precisely identified numerically in 
the present work. The plot of figure 19 confirms analytical results from [31] both 




Figure 19. Second cumulant of g as a function of first cumulant of g showing the 
universal behavior in the metallic regime. Different curves correspond to different 
transverse lengths (L y = 10,40,80) or magnetic disorder (J = 0,0.2). 

qualitatively in the shape of the curves and quantitatively in the values of fluctuations 
in both universality classes. In our study, values of UCF are reached with a maximal 
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Figure 20. Second cumulant of g as a function of longitudinal length, for L y = 80. 
Dots represent data for J — and squares for J = 0.2. The value of UCF is shown in 
each symmetry class (with or without magnetic disorder). UCF regime is reached in 
both cases. 




Figure 21. Second cumulant of g as a function of longitudinal length, for L y = 40. 
Dots represent data for J = and squares for J = 0.2. The value of UCF is shown 
in each symmetry class (with or without magnetic disorder. UCF regime is reached in 
the case of magnetic impurities but not for scalar impurities. 

error of 1% for GOE and 3% for GUE with respect to the analytical value of the UCF 
in the regime independent of (g) (i.e with much higher precision than e.g [38] and [30]). 

Let us finally comment the work of Z. Qiao et al. [33] who have performed a similar 
numerical Landuer study of ID transport in various universality classes, focusing mostly 
on the metallic regime. While both our study agree on the UCF (although we have a 
higher accuracy for /3 = 1), we did not find evidence for a second universal conductance 
plateau. This result would definitely deserves further study. 
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Figure 22. Variance of g as a function of longitudinal size. UCF are shown. Different 
curves correspond to different values of magnetic disorder J (J — 0.025 — > 0.4). 
Transverse length L y = 40. The only free parameter in analytical fits is the magnetic 
length L m . 
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Figure 23. Variance of g as a function of L x /L m . Transverse length L y = 40. 



5.2. The average conductance 

The main contribution to the average conductance is of classical origin. However a 
weak localization correction must be added when the quantum behavior of electrons is 
taken into account [12]. This quantum part manifests itself in the magneto-conductance 
behavior of long wires (larger than the the phase coherence length) where a weak 
magnetic field is sufficient to destroy this quantum correction by dephasing the various 
diffusing path with respect to each other [12]. We can write this average conductance as 
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Figure 24. Evolution of classical conductance g c iass = (g) — Sg with the longitudinal 
length for L y = 40 and J = 0.05, 0.1, 0.2, 0.4. Numerical data are well fitted using the 
conductivity as the only free parameter, different data correspond to different values 
of magnetic disorder. In inset is shown the set of numerical data for J — 0.2. The 
dotted line is the analytical fit without taking into account the contact conductance. 
Plain line is the same analytical fit as in the main plot. 



(without any magnetic field): 

xi Lfn) Qclass 

(J, L x ) + 5gwL(L x , L m ), (18) 

where g c i aS s is the classical part of the conductance and Sgwi is the weak localization 
correction. For a quasi one dimensional system the quantum correction reads the simple 
form [12] 
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The knowledge of the magnetic length L m we gained in the previous study of 
conductance fluctuations can now be used to completely characterize this weak 
localization contribution to the conductance. By subtracting the corresponding 
contribution to the average conductance, we obtain the classical conductance, plotted 
in figure 24 as a function of longitudinal length. Taking into account the contact 
resistance [39] in the two terminal setup, the expected expression for this classical 
conductivity reads 

9class(J, L x ) = 1 L -. (20) 

L v + a {J) 

Figure 24 shows this expression plotted for different values of magnetic disorder. The 
corresponding fitting parameter o"o is plotted as a function of J in figure 25. To compare 
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Figure 25. Evolution of conductivity oq with magnetic disorder J for L y = 40. Dots 
are numerical data. Plain line is the theory given by Einstein relation and Matthiesen 
rule for the conductivity. Agreement is good especially at low magnetic disorder. 



these results with theory, consider the Einstein relation which links the (Einstein) 
conductivity to the diffusion constant: 

cr = se 2 p (e F )D, (21) 

where s is the spin degeneracy and Po(^f) the electronic density of states at the Fermi 
level. By definition, the diffusion coefficient reads, for non magnetic impurities: 

D = v 2 F T e , (22) 

with Vp the Fermi velocity and r e the elastic scattering time. It can be related to the 
scalar disorder by: 

Te = - 1 2 (23) 

2irp riiVQ 

where rij is the impurity density and v 2 = W 2 /12. For more than one diffusive process, it 
is compulsory to use the Matthiesen rule that modifies scattering time r e in the following 
way: 

111 

-->- + — (24) 

7"e 7~e I'm 

where r m = L 2 m jD and is related to the magnetic disorder: 

1 



(25) 



Using this allows one to get the J dependance of the Einstein conductivity: 

a (J = 0) 

a 0(J) = i j 3_ J2 ' ( 26 ) 

In figure 25 we compare this expression with numerical evaluation of the conductivity. 
The good agreement between both curves provides an additional check of the 
correct determination of the magnetic dephasing length L m . Including the magnetic 
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Figure 26. Evolution of magnetic dephasing length L m with magnetic disorder J for 
L y = 40. Dots are numerical data extracted from the study of UCF and diamonds are 
numerical data extracted from the study of conductance correlations for two different 
spin configurations. Plain and dotted lines are analytical fits from perturbation theory 
at second order in J. Error bars are smaller than dots and diamonds sizes. 



disorder dependance of the diffusion coefficient (through Matthiesen rule), we obtain a 
perturbative expression to second order in J for this magnetic dephasing length: 

L m (J) = ^/D{J)r m {J) oc 1 (27) 

j w 2 , J 2 

J y -w + — 

On figure 26 we have plotted the numerical evaluation of the magnetic length as a 
function of magnetic disorder and the corresponding fit with eq.(27). We notice that it is 
also possible to obtain the magnetic length via the study of correlations of conductance, 
i.e via the study of (g(V, {S^ }i)g(V, {S!j- }i)) c , which goes beyond the scope of the 
present paper [11]. 

5.3. The third cumulant 

Finally we consider the third cumulant of the distribution of conductance. According 
to the analytical study of [31], this cumulant decays to zero in a universal way as (g) 
increases. Here in figure 28 we find a dependance of this decrease on the symmetry class: 
for GOE {g 3 ) c goes to zero in a monotonous way whereas it decreases, changes its sign 
and then goes to zero in GUE case. For (g) > 4 numerical errors are dominant, then this 
part of the curve is irrelevant. Moreover, for GUE this decrease seems to be universal 
whereas it depends on the transverse length for GOE. On figure 27, is represented 
the convergence of the skewness when increasing the number of configurations used to 
perform averages Nd es for both GOE and GUE. Plots show a good enough convergence 
of averages to conclude that the third cumulant of conductance is not zero for all values 
of (g). Notice that the maximal number of averages is 50000. Moreover this fast 
vanishing of the third cumulant confirms the faster convergence of the whole distribution 
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Figure 27. Plot of (g } c as a function of (g) in the metallic regime, averages are 
performed with various a number of configurations Nd es - Convergence curves are 
shown for L„ = 10 and J = or J = 0.2. 
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Figure 28. Plot of (g 3 ) c as a function of (g) in the metallic regime for L y = 10, 80 
and J = 0,0.2. 



towards the gaussian, compared to what happens in the insulating regime. Based on 
our numerical results, we cannot confirm nor refute the expected law (g 3 ) c oc l/{g) n , 
with n = 2 in GOE and n = 3 in GUE [40, 41]. 

6. Conclusion 

To conclude we have conducted extensive numerical studies of electronic transport in 
the presence of random frozen magnetic moments. Comparing and extending previous 
analytical and numerical studies, we have identified the insulating and metallic regimes 
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described by the universality classes GOE and GUE. We have paid special attention 
to the dependance on this symmetry of cumulants of the distribution of conductance 
in both metallic and insulating universal regimes. In particular, we have identified 
with high accuracy the domain of universal conductance fluctuations, and determined 
its extension in the present model. We have also determined precisely the so-called 
magnetic length L m which represents the elastic scattering length of the spin on magnetic 
impurities. This length is of primary importance in experiments as it controls the 
crossover between universality classes. This work paves the way for further studies 
of transport in metals with frozen magnetic impurities as we have clearly identified 
the range of the parameters to access the experimentally relevant metallic diffusive 
regime. One possible extension consists in considering evolution of the statistics of 
conductance as the magnetic disorder is varied, e.g. by rotating or flipping the spins of 
impurities. Comparing the conductance obtained in both spin configurations mimics 
the experimental measurement of the conductance of a low temperature canonical 
spin glass after two successive quenches [42, 11], without the necessary restrictions 
of analytical approaches [11]. Experimentally, this approach could give access to 
fundamental properties of a spin glass, that have never been measured. 

We thank X. Waintal for useful discussions. This work was supported by the 
ANR grants QuSpins and Mesoglass. All numerical calculations were performed on the 

computing facilities of the ENS-Lyon calculation center (PSMN). 

[1] Juliette Billy, Vincent Josse, Zhanchun Zuo, Alain Bernard, Ben Hambrccht, Pierre Lugan, David 
Clement, Laurent Sanchez-Palencia, Philippe Bouyer, and Alain Aspect. Direct observation of 
anderson localization of matter waves in a controlled disorder. Nature , 453:891-894, 2008. 

[2] Giacomo Roati, Chiara D'Errico, Leonardo Fallani, Marco Fattori, Chiara Fort, Matteo Zaccanti, 
Giovanni Modugno, Michele Modugno, and Massimo Inguscio. Anderson localization of a non- 
interacting bose-einstein condensate. Nature , 453:895-898, 2008. 

[3] P.G.N, de Vegvar, L.P. Levy, and T.A. Fulton. Conductance fluctuations of mesoscopic spin 
glasses. Phys. Rev. Lett. , 66:2380, 1991. 

[4] J. Jaroszynski, J. Wrobel, G. Karczewski, T. Wojtowicz, and T. Dietl. Magnetoconductance noise 
and irreversibilities in submicron wires of spin-glass n+- cdmntc. Phys. Rev. Lett. , 80:5635, 
1998. 

[5] G. Neuttiens, J. Eom, C. Strunk, H. Pattyn, C. Van Haesendonck, Y. Bruynseraede, and 
V. Chandrasekhar. Thermoelectric effects in mesoscopic aufe spin-glass wire. Europhys. Lett. , 
42:185, 1998. 

[6] T. Capron, A. Perrat-Mabilon, C. Peaucelle, T. Mcunicr, D. Carpentier, L. P. Levy, C. Baeuerle, 
and L. Saminadayar. Remanence effects in the electrical resistivity of spin glasses. EPL , 93(2), 
JAN 2011. 

[7] E. Abrahams, P.W. Anderson, D.C. Licciardello, and TV. Ramakrishnan. Scaling theory of 
localization: Absence of quantum diffusion in two dimensions. Phys. Rev. Lett. , 42:673, 1979. 

[8] F. Evers and A. D. Mirlin. Anderson transitions. Rev. Mod. Phys. , 80:1355, 2008. 

[9] Shinsei Ryu, Andreas P Schnydcr, Akira Furusaki, and Andreas W W Ludwig. Topological 
insulators and superconductors: tenfold way and dimensional hierarchy. New Journal of Physics , 
12(065010), 2010. 

[10] A. D. Mirlin. Statistics of enrgy levels and eigenfunctions in disordered systems. Phys. Rep. , 
326:259, 2000. 

[11] G. Paulin, T. Capron, D. Carpentier, T. Mcunicr, C. Bauerle, L.P. Levy, and L. Saminadayar. 
Spin flip induced conductance fluctuations, preprint, 2011. 



Cross-Over between universality classes in a magnetically disordered metallic wire. 26 

[12] E. Akkcrmans and G. Montambaux. Mesoscopic Physics of electrons and photons. Cambridge 
University Press, 2007. 

[13] M. Mezard, G. Parisi, and M. Virasoro. Spin Glass Theory and Beyond . World Scientific, 1987. 
[14] O. N. Dorokhov. Zh. Eksp. Teor. Fiz. , 85:1040, 1983. 

[15] P. A. Mello, P. Pereyra, and N. Kumar. Macroscopic approach to multichannel disordered 

conductors. Ann. Phys. , 181:290-317, 1988. 
[16] R. Landauer. Spatial variation of currents and fields due to localized scatterers in metallic 

conduction. IBM J. Res. Dev. , 1:223, 1957. 
[17] L. I. Deych, M. V. Erementchouk, and A. A. Lisyansky. Scaling in the one-dimensional anderson 

localization problem in the region of fluctuation states. Phys. Rev. Lett. , 90:126601, 2003. 
[18] M. Buttiker, Y. Imry, R. Landauer, and S. Pinhas. Generalized many-channel conductance formula 

with application to small rings. Phys. Rev. B, 31:6207-6215, 1985. 
[19] D. S. Fisher and P.A. Lee. Relation between conductivity and transmission matrix. Phys. Rev. 

B, 23:6851-6854, 1981. 

[20] A. Croy, R. A. Roemer, and M. Schreiber. Localization of electronic states in amorphous materials: 
recursive green's function method and the metal-insulator transition at ej^O. In K. Hoffman and 
A. Meyer, editors, Lecture notes in Computational Science and Engineering , volume Lecture 
notes in Computational Science and Engineering. Springer (Berlin), 2006. 

[21] C. W. J. Becnakker. Random-matrix theory of quantum transport. Rev. Mod. Phys. , 69:731, 
1997. 

[22] K. Slevin, P. Markos, and T. Ohtsuki. Reconciling conductance fluctuations and the scaling theory 

of localization. Phys. Rev. Lett. , 86(16):3594, 2001. 
[23] P. Markos. Weak disorder expansion of lyapunov exponents of products of random matrices: A 

degenerate theory. Journal of Statistical Physics , 70(3/4):899, 1993. 
[24] J-L. Pichard. Random matrix theory of scattering in chaotic and disordered media. In P. Sebbah, 

editor, Waves and Imaging through Complex Media . Kluwer Academin Publishers, 2001. 
[25] J-L. Pichard In NATO ASI SERIES B254 B. KRAMER, editor, Quantum Coherence in 

Mesoscopic systems , page 369. PLENUM, NEW YORK, 1991. 
[26] I. Larkin. Zh. Eksp. Teor. Fiz. , 85:764, 1983. 

[27] I.V. Lerner and Y. Imry. Magnetic-field dependence of the localization length in anderson 

insulators. Europhys. Lett. , 29:49-54, 1995. 
[28] A. MacKinnon and B. Kramer. One-parameter scaling of localization length and conductance in 

disordered systems. Phys. Rev. Lett. , 47(21):1546, 1981. 
[29] K. A. Muttalib and P. Wdlfle. "one-sided" log- normal distribution of conductances for a disordered 

quantum wire. Phys. Rev. Lett. , 83(15):3013, 1999. 
[30] P. Markos. Dimension dependance of the conductance distribution in the nonmetallic regime. 

Phys. Rev. B , 65:104207, 2002. 
[31] L.S. Froufc-Perez, P. Garcia-Mochales, P.A. Serena, P.A. Mello, and J.J. Saenz. Conductance 

distributions in quasi-one-dimensional disordered wires. Phys. Rev. Lett. , 89(24):246403, 2002. 
[32] K. A. Muttalib, P. Wolfle, A. Garcia-Martm, and V. A. Gopar. Nonanalyticity in the distribution 

of conductances in quasi-one-dimcnsional wires. Europhys. Lett., 61:95, 2003. 
[33] Z. Qiao, Y. Xing, and J. Wang. New universal conductance fluctuation of mesoscopic systems in 

the crossover regime from metal to insulator. ArXiv:0910.3475vl , 2009. 
[34] V. A. Gopar, K. A. Muttalib, and P. Wdlfle. Conductance distribution in disordered quantum 

wires: Crossover between the metallic and insulating regimes. Phys. Rev. B , 66:174204, 2002. 
[35] C. W. J. Beenakker and H. van Houten. Quantum transport in semiconductor nanostructures. 

Solid State Physics , 44:1, 1991. 
[36] M.R. Zirnbauer. Super fourier analysis and localization in disordered wires. Phys. Rev. Lett. , 

69:1584, 1992. 

[37] A. D. Mirlin, A. Mueller-Groeling, and M. R. Zirnbauer. Conductance fluctuations of disordered 
wires: Fourier analysis on supersymmetric spaces. Ann. Phys. , 236:325, 1994. 



Cross-Over between universality classes in a magnetically disordered metallic wire. 27 



[38] M. Cieplak, B.R. Bulka, and T. Dietl. Universal conductance fluctuations in spin glasses. Phys. 
Rev. B , 44:12337, 1991. 

[39] H. L. Engquist and P.W. Anderson. Definition and measurement of the electrical and thermal 

resistances. Phys. Rev. B , 24:1151, 1981. 
[40] A. M. S. Macedo. Random-matrix approach to the quantum-transport theory of disordered 

conductors. Phys. Rev. B , 49:1858, 1994. 
[41] M.C.W. van Rossum, Igor V. Lerner, Boris L. Altshuler, and Th. M. Nieuwenhuizcn. Deviations 

from the gaussian distribution of mesoscopic conductance fluctuations. Phys. Rev. B , 55:4710, 

1997. 

[42] D. Carpentier and E. Orignac. Measuring overlaps in mesoscopic spin glasses via conductance 
fluctuations. Phys. Rev. Lett. , 100:057207, 2008. 



